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Abstract 

We compute the nucleon matrix elements of QCD operators of twist 3 and 4 in the 
instanton vacuum. We consider the operators determining 1/Q^-power corrections 
to the Bjorken, Ellis-Jaffe and Gross-Llewellyn-Smith sum rules. The nucleon is 
described as a soliton of the effective chiral theory derived from instantons in the 
l/A'^c~expansion. QCD operators involving the gluon field are systematically rep- 
resented by effective operators in the effective chiral theory. We find that twist-3 
matrix elements are suppressed relative to twist-4 by a power of the packing fraction 
of the instanton medium. Numerical results for the spin-dependent {S'^\ f^"^^) and 
spin-independent twist-3 and 4 matrix elements are compared with results of other 
approaches and with experimental estimates of power corrections. The methods 
developed can be used to evaluate a wide range of matrix elements relevant to DIS. 
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1 Introduction 



The discovery of scaling in deep-inelasting scattering experiments can be seen as the 
starting point of modern strong interaction physics, giving rise to the parton model and 
supporting the idea of asymptotic freedom. QCD was established as a theory which 
describes logarithmic corrections to scaling in the asymptotic region. The predictions of 
perturbative QCD for the scale dependence of structure functions in leading and next- 
to-leading order have been well confirmed by experiment However, in general the Q^- 
dependence of structure functions is subject to power corrections, which lie outside the 
scope of perturbative QCD. These corrections are governed by scale parameters associated 
with non-perturbative vacuum fluctuations. 

The standard tool to analyze the scale dependence of structure functions in QCD has 
been the operator product expansion. The part of nucleon structure functions exhibiting 
only logarithmic scale dependence is determined by nucleon matrix elements of operators 
of leading twist, while a class of power corrections are associated with operators of non- 
leading twist [0]. In a partonic language twist-2 operators are one-particle operators 
counting the contribution of individual partons to the structure functions, while higher- 
twist operators can be interpreted as measuring the effect of the chromoelectric and - 
magnetic field of the nucleon, or effects of correlation between partons, on the structure 
functions. For a quantitative description of power corrections one needs the values of the 
matrix elements of these operators. These can be regarded as fundamental characteristics 
of the nucleon. Clearly, an estimation of these quantities requires a theory of the non- 
perturbative gluon degrees of freedom. 

A microscopic picture of the non-perturbative fluctuations of the gluon field is provided 
by the instanton vacuum. The relevance of instantons to the hadronic world was explored 
already before a quantitative theory of the instanton medium existed [|], |], ^ . A treatment 
of the interacting instanton ensemble by means of the Feynman variational principle 
showed how the instanton medium stabilizes itself . The most prominent characteristic 
of the instanton vacuum is the small packing fraction of the medium, i. e., the small ratio of 
the average size of the instantons in the medium to the average separation between nearest 
neighbors. A value p/R ^ 1/3 was found both from phenomenological considerations 
0] and from the variational principle 0. This small parameter is the starting point 
for a systematic analysis of non-perturbative phenomena in this picture, such as chiral 
symmetry breaking, vacuum condensates etc. We note that the instanton vacuum has also 
been studied in lattice experiments using the cooling method. The measured properties of 
the instanton medium found after cooling of the quantum fluctuations are close to those 
computed from the variational principle [^. 

The main success of the instanton vacuum is its explanation of the dynamical break- 
ing of chiral symmetry, which is the most important non-perturbative phenomenon de- 
termining the structure of light hadrons, including the nucleon. The mechanism of chiral 
symmetry breaking is the delocalization of the fermionic zero modes associated with the 
individual instantons in the medium |0, |^, resulting in a finite fermion spectral density 
at zero eigenvalue, which is proportional to the chiral condensate |^. In an alternative 
formulation of chiral symmetry breaking one derives the effective action of fermions in the 
instanton medium by integrating over the instanton coordinates in the ensemble. It has 
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the form of a Nambu-Jona-Lasinio model |Ty] with a many-fermionic interaction with a 
specific spin-flavor structure, the so-called 't Hooft determinant 0. Due to this interac- 
tion the quarks develop a dynamical mass, and a Goldstone pion appears as a collective 
excitation. The effective theory of massive quarks and pions applies to quark momenta up 
the inverse instanton size, ~ 600 MeV, which acts as a natural cutoff. The ratio of the 
dynamically generated quark mass to the cutoff is directly related to the packing fraction 
of the instanton medium, Mp oc (p/R)'^. Thus, the diluteness of the instanton medium 
ensures that the picture of interacting "constituent" quarks applies in a parametrically 
wide rangeQ. 

An immediate application of the effective chiral theory derived from instantons is the 
quark-soliton model of baryons |T^. In the limit of a large number of colors, Nc, the 



nucleon can be viewed as Nc "valence" quarks bound by a self-consistent hedgehog-like 
pion field. This approach gives a very reasonable description of a variety of nucleon and 
A properties such as the A^-A splitting, electric formfactors, magnetic moments, axial 
coupling constants etc. |T^. Recently it was shown that this approach also allows to 



compute the twist-2 quark distribution functions of the nucleon at a low normalization 
point of order p~^ |15|, |16|, O. 



The instanton vacuum, with the resulting effective chiral theory, allows to evaluate 
hadronic matrix elements of QCD operators involving the gluon field. In ref.[]T8[ a method 
was developed by which "gluonic" operators can systematically be represented as effective 
many-fermionic operators in the effective chiral theory, the hadronic matrix elements of 
which can then be computed using standard techniques. This "fermionization" of QCD 
operators is possible relying entirely on the approximations which are already inherent 
in the effective theory, namely the diluteness of the instanton medium and the l/Nf.- 
expansion — no additional assumptions are required. It was shown in |TB| that this 
approach preserves the essential renormalization properties of QCD; for example, the 
QCD trace and U{1) anomalies are realized in this approach at the level of hadronic 
matrix elements. This method is thus well suited for computing the matrix elements of 
the QCD operators which arise in the description of deep-inelastic scattering, both for 
leading and non- leading twist. 

In this paper we evaluate the nucleon matrix elements of the operators of twist 3 and 
4, which describe power corrections to the lowest moments of polarized and unpolarized 
nucleon structure functions, using the approach of [0. Our aims here are manifold. 



First, we want to develop the methods for computing nucleon matrix of operators of 
the type arising in the description of power corrections to DIS, and demonstrate the 
consistency of this approach; for example, we verify that the QCD equations of motion 
are preserved by the effective operators. Second, we want to classify the various higher- 
twist matrix elements according to their magnitude with regard to the parameters of our 
approach, the instanton packing fraction, Mp oc [p/R)"^, and the formal parameter l/Nc- 
We shall see that a characteristic difference between operators of highest twist (lowest 
spin) and lower twists (higher spins) emerges as a consequence of the 0(4)-symmetry 
of the instantons. We also investigate the role of the non-perturbative gluon degrees of 
freedom in twist-2 operators within our approach. Finally, we shall obtain numerical 



^For a recent review of the instanton vacuum see |T^, [l2| . 
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values for the higher-twist matrix elements and compare with experiment as well as with 
other theoretical estimates. The power corrections are sensitive to the form of the non- 
perturbative vacuum fluctuations, so their study and comparison with experiment will 
shed some light on the srtructure of these fluctuations. 

The plan of this paper is as follows. In section 2 we discuss the twist-3 and 4 operators 
whose matrix elements determine power corrections to the lowest moments of polarized 
and unpolarized nucleon structure functions. In section 3, we first summarize the proper- 
ties of the instanton vacuum and the effective chiral theory. We then derive the effective 
operators for the QCD operators of leading and non-leading twist which generally arise 
in the description of nucleon structure functions. In the next step we compute the matrix 
elements of the twist-3 and 4 operators in quark states and discuss the parametric sup- 
pression of twist 3 relative to twist 4 in the instanton packing fraction. Furthermore, we 
show explicitly that the QCD equations of motion are preserved by the effective operators. 
In section 4 we briefly describe the nucleon in the effective chiral theory. We then derive 
the expressions for the matrix elements of the gluonic operators of interest. We clas- 
sify the spin-dependent and independent higher-twist matrix elements according to their 
magnitude in the l/A^^c-expansion. We then compute the A'^c-leading higher-twist matrix 
elements and give crude estimates for the subleading ones. A discussion of the numerical 
results and comparison with experimental numbers and results of other non-perturbative 
methods are presented in section 5. Section 6 offers conclusions and an outlook. 



2 Power corrections to DIS 



The Q^-dependence of nucleon structure functions in QCD is usually investigated using 
the operator product expansion Expanding the product of electromagnetic or weak 
currents in the forward Compton amplitude at near light-cone separations in operators 
of increasing twist, one can identify the scaling and power-suppressed parts of moments 
of the structure functions. The leading (l/Q^-) power corrections to the isovector and 
isosinglet combinations of the first moment of the polarized structure function gi (the 
Bjorken |jT9| and EUis-Jaffe pO| sum rules) are given by^ 



/ dx gi~^^{x,Q^ 
Jo 



+ 
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Ml 

2 («. 



{a% + 4d' 



(2) 
NS 



+ 
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9Q 



(2) 
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+ Adf ' ^ 



(2.1) 
(2.2) 



Here M^r denotes the nucleon mass. The constants a^^\a^'^^ are defined by the matrix 
elements of the twist -2 operators of spin 1 and 3, respectively. More generally. 



{PSn\iPT\^j,Vp, . . . V(3„}i:\PSn) - traces 

= -2Mn{2Ts) a^^l S{a • • • P/3n} - traces, 
(PS* I V^7{c, 75 V/3i . . . V/3„}^|PS') - traces 



(2.3) 



^ We note that there is a mistake in ref. concerning the coefficients in this formula, see ^ 



5 



= -2Mn S{aPi3i ■ ■ ■ Ppn} - traces, 
(n = 0,2,...). 

The braces denote complete symmetrization in the respective indices, 
nucleon 4-momentum, 5* the nucleon polarization vector, 



(2.4) 



Here P is the 



0, 



-1. 



(2.5) 



The minus sign in eqs.( p.3| , p.4| ) relative to ref.|^2[ comes because our 75 differs from the 
one of Bjorken and Drell by a minus sign; this definition is used in refs.[^ |T3|, [T^ and also 
ref.[^. In particular, and a^"* are identical to the nucleon isovector and isosinglet 
axial coupling constants 



''NS 



(3) 

9a , 



(0) 

9 a ■ 



(2.6) 



In eg. (p.3[) T3 = ±1/2 for proton and neutron, respectively, and denotes the Pauli 
matrix in flavor indices. In the following it will always be understood that the flavor- 
nonsinglet (NS) and singlet {S) parts of bilinear operators ijj . . .ip and their matrix ele- 
ments are defined as in eqs. (|2.3| , |2.4| ); we shall not write the flavor dependence of matrix 
elements explicitly. (We consider here the case of two light flavors; the generalization of 
the formulas here and below to S't/(3)-flavor is straightforward.) 

Of the power corrections in eqs. (|2.1| , |2.2|) , the terms proportional to the twist-2 matrix 
element a*-^-* originate from the expansion of Nachtmann moments (target mass correc- 
tions) im, while the other terms represent contributions of operators of non- leading twist. 
The constant c?'-^^ parametrizes the matrix element of the twist-3 spin-2 operator 



{PS\tP (7°F^^ + 7^F"^) ^\PS)- traces 
= 2Mn d'^^^ hP^'P'^S^ - P^pf^S"" - P'^P^S'^ + {a^ p)- traces 



(2.7) 



for both flavor non-singlet and singlet, see eqs.( p^.3| , O). It is understood that trace terms 
are to be subtracted on both sides, since one is interested only in the twist-3 part. The 
same operator contributes also to the non-power suppressed-part of the third moment of 
the structure function (72, 



Jo 



, 1 >(2) ,r,(MR 



1 

--a 



(2) 



Finally, Z^^-* is given by the matrix element of the twist-4 spin-1 operators p2 

{PS\^PjJ''''^\PS) = 2M3/(2)5^ 
again for both flavor non- singlet and singlet. 



(2.8) 
(2.9) 

(2.10) 
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In unpolarized neutrino scattering there are l/Q^-power corrections to the first mo- 
ment of the isovector parity-even structure function Fi and the isosinglet parity-odd 
structure function F3 (Gross-Llewellyn-Smith sum rule 




3Q2 



8M 



9Q 



NJ2) ^ 
NS ^ ^ 



1 - 



27Q 



2^5 



(2.11) 
(2.12) 



Here, C(2) is determined by the spin-independent nucleon matrix element of the twist-4 
spin-1 operator^ [p4| , 



(2.13) 



spin 



We note that in both the polarized and unpolarized case, eqs.(P?T|, |2.2D and eqs.( |2.11 
|2.12|) , there are, of course, logarithmic corrections, which we have not written down |[l[]. 
They originate from the scale dependence of the coefficients of the OPE, which is described 
by the renormalization group equation. 

In this paper we evaluate the nucleon matrix elements of the twist-3 and 4 operators, 
eqs.( p.7| , p.lO| ) and eq.( |2.13| ), in the instanton vacuum. The matrix elements depend on 



the normalization point of the operators, which in our approach is of the order of the 
inverse average instanton size, p^^ ~ 600 MeV. Values at other normalization points can 
easily be obtained using the renormalization group equation. 

We emphasize that we are not dealing here with the instanton-induced contributions 
to DIS which arise from non-perturbative power corrections to the coefficient functions of 
the OPE. Those are the leading semiclassical corrections associated with small instantons 
of size p ~ 1/Q [BR 



3 From QCD to the effective chiral theory 

3.1 Effective chiral theory from the instanton vacuum 

The basis of our description is the medium of independent instantons with an effective size 
distribution, which was obtained by Diakonov and Petrov as a variational approximation 
to the interacting instanton partition function ||5[. In the large- Ai'c limit the width of 
the effective size distribution is of order 1/6 ~ 1/^c, so we may assume instantons and 
antiinstantons of fixed size, 

p ~ (600MeV)-\ (3.1) 

It is sufficient to consider a medium with equal number of instantons and antiinstantons, 
= = N/2, since we shall not be concerned with topological fluctuations or effects 
related to the f/(l)-anomaly here. 

^We follow here ref.Hl, with ((O^'^''-^)) = M'^c^g^j^g. Note that the conventions of differ from 
those of ; see the corresponding remarks in ref . |Q . 
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Fermions couple to the instantons and antiinstantons (J's and J's in the following) 
mainly through the zero modes. In the zero-mode approximation, which is justified by 
the diluteness of the instanton medium, the interaction of the fermion field with given 
fiavor / with one / (/), with collective coordinates z (center) and lA (color orientation) is 
given byfl [0, |TB 



(27r)4 7 (27r)' 

1 1±75 

-7k7a — - — 



(3.2) 



(We use the standard ip"^ = iip for the quark field.) Here F{k) is a form factor proportional 
to the Fourier transform of the wave function of the fermion zero mode, which drops to 
zero for momenta of order p~^, 



m 



-tj^[I,{t)Ko{t) - h{t)K^{t)] 



t 



t 

oo 



(3.3) 



Ik p. 



The matrices are Nc x Nc matrices with (r, =Fi) in the upper left corner and zero 
elsewhere. 

The effective fermion action {Nf flavors), which is derived by integrating over the 
coordinates of the instantons in the medium and flnding the saddle point of the fermion 
integral in the large-iVc limit [^O], 0, is of the form 

Here, Y± denote the 2iVj-fermion vertices 

f N M\^' r r 



(3.4) 



47r2p2y 



iM)^f J d^x det J±(a;) 



(3.5) 



where J±{z)fg are color singlet currents which are Nf x A'^^-matrices in flavor, 

d^k 



J±{x) 



fg 



{2tiY J (27r) 



^'\e-^('=-')-F(fc)F(/) ^}(A:)i4^^,(/), 



(3.6) 



and the determinant in eq.( |3.5| ) is over flavor indices. The vertices eq.( p.5| ) have the 
spin-flavor structure of the 't Hooft determinant [Q]. The dynamical quark mass, M, 
is determined by the saddle-point equation. Its square is proportional to the instanton 

^To avoid confusion with the pion field we denote in this paper the color orientation matrices of the 
instanton in the fundamental representation of SU{Nc) by not by U as in ref . p8| . 
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density, since chiral symmetry breaking is a collective effect involving all instantons in the 
medium. Thus, parametrically. 



(Mp) 



OC 



OC — 



(3.7) 



The diluteness of the instanton medium thus ensures that the dynamical quark mass is 
small compared to the momentum cutoff, p~^, so that the effective theory of massive 
quarks described by eq.( p.4|) is valid in a parametrically wide range. This fact is central 
to the concept of the effective chiral action. 

In the case of one quark flavor, Nf = 1, the 't Hooft vertex, eq. (|3.5|) , is simply a mass 
term, and the effective fermion action becomes 



'S'eff[V'^^]JV/= 



^\k) ^ -iMF\k) ilj{k). 



(3.8) 



For more than one quark flavors, Nf > 1, the effective action, eq.( t3.4|) , describes a system 
of quarks with many-fermionic interactions. The partition function of this theory can be 
written in bosonized form by introducing meson fields [jlO|, |T8 |. 



Vtt I Vilj^Vip exp / (fx 



'Nf 



Nf 



+iMY.[U+{x)f,J^ 

f,9 



fg 



+ U-{x)fgJ4x 



}f9\ 



(3.9) 



Here, 7r(x) is the pion field, and 

U+{x) = U{x), 
U{x) = exp [z7r'^(x)r° 



U_{x) = U^{x), 



(3.10) 



We have retained in eq. ( p.9|) only the Goldstone degrees of freedom; other mesons {e.g. 
rho, sigma) have masses of order of the cutoff, p~^, and thus should not be considered 
as dynamical degrees of freedom in this effective theory. Eq.( [3.9|) describes the minimal 
chirally invariant coupling of quarks to Goldstone bosons. 

Eq.( p.9| ) is the starting point for the description of the nucleon, which in the large- A'^c 
limit is described by a classical pion field of hedgehog form, see Section |. 



3.2 Effective gluon operators 

The higher-twist operators arising in the description of power corrections to DIS explicitly 
involve the gluon field. Following the approach of [|TB| we now determine the effective 
fermion operators representing these operators in the effective chiral theory. 

In a dilute medium of instantons the interaction of a gluonic operator with the light 
fermions is dominated by the contribution of a single /(/); many-instanton contributions 
are suppressed by powers of the packing fraction. Thus, when constructing the effective 
operator corresponding to a given function of the gauge field, one evaluates this function 
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in the field of one /(/)• Tlie functions of the gauge field appearing in the QCD operators 
of twist 2, eqs. (^l3| , p.4|) , and twist 3 and 4, eqs. (|2.7| , |2.10| , |2.13| ), have the property that 
in the field of one /(/) they are proportional to the 't Hooft eta symbol. For instance, 
the field strength and dual field strength of one I (/) in singular gauge have the form 

n2 



^l|3 



where (r/-);^^ = f]"^^, (r/+);^^ 
given by 



±i^|,(^)/(/), 

rj'^j^. The same applies to the gauge field itself, which is 



(3.11) 
(3.12) 



2p^ 



(x^ + p^)x' 



(3.13) 



or, more generally, to the products of gauge fields which arise when multiplying out the 



product of covariant derivatives in the twist-2 operators of arbitrary spin, eqs. (|2.3| , |2.4| ), 

~2 



{Alix) . . . AXix) 



tr 



(')^)l/±,^.«,...»„(a:), (3.14) 



where f±^v/3^...i3^{x) denotes a Lorentz tensor depending on x whose precise form can be 
determined by explicitly computing the product using eq.( p.l3| ). 

To account for all cases of interest^, we consider in the following a general "gluonic" 
operator of the form 



A" 

Oai...a,/3i.../3,(a;) = V^l^^) yF^i ...a, ^ (x) J^^^ ...^^ (x) , 



(3.15) 



[X 



'F'^^...pJ[A{x)] denotes a 



where Ta^^ ^a,- is a matrix in Dirac spinor indices and f. 
function of the gauge field in the adjoint representation, which in the field of one /(/), in 
standard orientation and centered at zero, is given by 

^l...pMm = iv^)%^±,,up....Prix), (3.16) 

with J-'±,fMui3i...i3r{x) are general tensor functions of the coordinate, x. For a general instan- 
ton with center z and color orientation given by the SU{Nc) matrix U, cf. eq. (|3.2|) , the 
function eq.( p.l6| ) takes the form 



^(3i.../3ri^)l{I) 



(3.17) 



Following [|T8| we define a 2iVj-fermion vertex as the average of the product of the function 
eq. (|3.17|) and the instanton-fermion interaction potential, eq. (|3.2| ), over the collective 
coordinates of one instanton, 



(>v±);,...;3J^^v^](^) 



N ) 



I 



Nf ^ Nf 

^d'zjdu :Fi...pM)iii) n y±m^^fi (3-18) 
/=i 

^We note that there are functions of the gauge field in the adjoint representation which in one /(/) 
are not proportional to the eta symbol. They are constructed with the help of the totally symmetric 
structure constants, d"'"^. These functions lead to effective operators of a form different from the one 
described below. We do not need to consider them here. 
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Performing the integral over color orientations in the leading order of l/Ai'c, as in eq.( |3.5| ), 
one obtains 



2V 



Nf~l 



j d*zT±^^^p,,„pXx- z) J2 ^±,/.^(^)/9 det'^gJ±(z). (3.19) 



/,9=1 



Here J±{z) are the left- and right-handed color singlet currents, eq.(|3.6|), while J±^^i, is 
the color-octet Lorentz tensor current, 



dH 



A" 



(27r)4 J (27r) 



F{k)F{l) ipUk) — 



a 



1±75 



fiu 



(3.20) 



where a^u = In eq.( |3.19[ ) detj^ denotes the minor of the determinant in which 

the /-th row and g-th column are omitted. The form of the vertex eq.( p.l9 ) is intuitively 
plausible: the "effective gluonic field" is given by the color-octet current of the quark 
fields convoluted with the function of the field of one instanton. Note that the vertex 
eq. (|3.19|) is chirally invariant — the tensor currents, eq. (|3.20|) , have the same chirality 
properties as the scalar-pseudoscalar ones, eq. (pl6|) , and with respect to flavor the terms 
in eq.( p.l9 ) combine to give the determinant, which is an invariant. 

The effective operator which represents eq.( |3.15| ) in the effective chiral theory is given 
by MM 



"O" 



Oil. ..a 



./3l.../3J^^^](a;) 



^^^^^^)_T^,,„^Mx) [(lV+)^,...^,(x)iv+7^+ + (lv_)^^...^^(x)iv_7^_J . (3.21) 



Here = = N/2 are the numbers of I's and I's in the ensemble. The 7l± are 
many-fermionic vertices which act as the "inverse" of the 't Hooft vertex, eq.( p.5| ), in 
saddle-point approximation (large- A^c limit); symbolically 



7^^ 



(3.22) 



(An exphcit representation for TZj-lip^ip] can be found in [@.) Eq. (|3.22| ) is to be under- 
stood in the following sense. The vacuum expectation value of the vertices is 



(3.23) 



while the connected average of TZ±. with other fermion fields is equivalent to the insertion 
of a 't Hooft vertex, eq.(p.5|). 



\ I conn 



(3.24) 



We note that the operator character of the vertices 7^±, eq.( p.22D , is relevant only when 
computing connected averages of fermion fields with effective operators which have a 
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non-zero vacuum expectation value, such as the operator ip^ip., see section p.4| . For 
operators with zero vacuum expectation value, which we shall mostly be concerned with 
in the following, the vertices TZ± in eq. ( ^.211 ) can be replaced by their vacuum expectation 
values, eq.( |3.23|) , so that the effective operator eq. ( |3.21|) reduces to 



"O' 



OL\...OL 



-A.../3JV'^v^](a;) 



(3.25) 



(Diagrams in which the factors TZ± are connected to external quark lines are suppressed 
by 1/Nc relative to the disconnected terms.) 

The effective operator eq. (|3.25| ) is given by the sum of / and / contributions. If 
the function of the gluon field, eq.( |3.16| ), takes the same value for / and /, J-'-,^ui3i...(3s = 
^+,fiui3i...i3s: as for example in the case of the field strength, F^,^, eq. (|3.11| ) or the gauge field. 



A^, eq.( |3.13D , the sum of vertices (Fjc-_|_)^^ is an operator of natural parity (true tensor). 

If ^±,iiuf3i...i3s takes opposite values for J's and J's, as for example for the dual field, F^,^, 
eq.( p.l2| ), the sum of vertices represents an operator of unnatural parity (pseudotensor) . 

In the case of one quark flavor, Nf = 1, the fermion vertex, eq. (|3.18|) , is bilinear in 
the fermion field, 



iM 



/ d^zJ^±,^,p,...p^{x - z)Jl^,{z). (3.26) 



where J±^y are the left- and right-handed currents, eq. ( |3.20| ). The vertex can be expressed 
in momentum representation as 



(>v±);,...;3Jv^^^](^)iv,=l 

_ iM f d'^ki f d^k 



27r 



a 



1±75 



(3.27) 



Here the Fourier transform of the function of the instanton field, eq. (|3.16| ), is defined as 



±,fj.iyf3i...f3, 



(k) 



J d^xe '^"^ J^±,^,u|3^...|3Sx). 



(3.28) 



The effective operator, eq. (|3.25|) , in this case is a four-fermionic operator. It has the 
form of a color-octet current-current interaction, with one current formed by the fermion 
fields present in the original QCD operator, the other by the non-local instanton-induced 
vertices. Graphically it can be represented as in Fig.([^). 

In the case of more than one quark flavor, Nf > 1, the effective operator is obtained at 
first as a (2 A^j + 2)-fermionic operator, eq. ( [3.251 ). When passing to the bosonized effective 
theory, eq. ( p.9|) , one may convert also the effective operator to a simple "bosonized" form. 
Consider a correlation function with the effective many-fermionic operator, eq. (|3.25 ), in 
the bosonized theory. The correlation function is evaluated by averaging first over the 
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fermion fields in the background of the pion field, later over the pion field (in saddle- 
point approximation). Let us denote the result of the averaging over the fermions fields 
by (. . .)7r. In the leading order of l/N^. one finds 

(. . . ^n^)r.,...„.y^(x) [{Y^rTp,...pM) + {Y.^Tp,...pS^) 



+ 



(3.29) 



since the currents J± appearing in the subdeterminant are color-singlet. (We remind 
that the J±, eq.( p.6|) , are matrices in flavor indices.) The ellipses here stand for other 
functions of the fermion fields, e.g. hadronic currents. Computing the last factor in 
eq.( p.29| ) explicitly one finds in leading order in Mp: 



Nf~l 



Nf-l 



U+(z 



(3.30) 



where U±{z) is the pion field, eg. ( P.10|) . The first equality holds in leading order of 
l/Nc, and in the last step we have used the self-consistency condition defining the saddle 
point of the effective chiral theory. Thus, in the bosonized effective theory, eq.( |3.9| ), 
the instanton-induced 2A^j-fermionic vertex, eq. (|3.18 ) is equivalent to the 2-fermionic 
"bosonized" vertex 



bosonized 



iM 



f,9 



(3.31) 



which is manifestly chirally invariant. This expression formally includes the case Nf = 1, 
eq. (|3.26|) , if the pion field is set to zero, U± = 1. Accordingly, the effective operator 
representing the gluonic operator, eq.( p.l5|) , is again equivalent to a 4-fermionic operator, 

"0"„,...«,;3,.../3J^^^;v^](x) 

= -#^(x)yr„^...„^V(x) [(y^+))^i.../3^(a;)bosonized + (>>-))^i.../3^ (a;)bosonized] , (3.32) 

which now, however, is coupled to the pion field. This representation of the effective 
operator is useful for calculations of nucleon matrix elements, see Section It is also of 
principal interest, since it makes manifest the chiral invariance of the effective operator. 



3.3 Matrix elements of higher— twist operators in quark states 

Before computing the nucleon matrix elements of higher-twist operators it is instructive 
to analyze their matrix elements in "constituent quark" states, i.e., the massive quark 



13 



states corresponding to the fermion fields of tlie effective chiral theory, eq.( p.4| ). A pole 
appears in the quark propagator after continuation to Minkowskian momenta, that is, at 
Euchdean = -M^ + 0{Mp). 

Let us first compute the quark matrix element of the effective operator corresponding 
to the general gluonic operator, eq.(PTTH). It is sufficient to consider one quark flavor. 



Nf = 1, in which case the effective operator is given by the 4-fermionic operator, eq. (|3.25| ) 



with the instanton-induced vertices eq.(^]^). It is convenient to carry out the entire 
calculation using Euclidean fields and Euclidean vector components, and to continue to 
= — at the end. The quark propagator in the effective theory is given by 

G(k) - (3 33) 

When computing the matrix elements of the 4-fermionic operator between quark states 
there are two contractions (see Fig.|l|). One obtains 



j F{p)F{p - k) 
\{p- kf + W^F^p - k) 

F{p)F{p + k) 
^ {p + ky + M^F\p + k) 



Ap,s r«,...„. U-lit+ iMF\p - k)) a, ^ + 



2 



+ — — 2~^) ■ (3.34) 

It is understood that one sets p^ = — after evaluating the integrals. Here 

A,,s = u{p,S)u{p,S) = (l + ^75>g) (3.35) 

is the projector on quark states with definite momentum and polarization vector, where 
p and S are Euclidean 4-vectors satisfying 

p2 = ^2 = 1, S-p = 0. (3.36) 

A factor of Nc resulting from the quark loop has canceled the factor 1/Nc in the effective 
operator, eq.( p.26| ), so that the result is 0{N^). Eq.( |3.34D is understood as the diagonal 



matrix element for a quark state of given color. We remind that formula eq. (|3.34|) is valid 
for operators with zero vacuum expectation value, see the discussion in Section \i.2\ 
We now apply the general formula eq. (|3.34| ) to the spin-dependent quark matrix el- 



ements of the twist -3 and 4 operators, eqs.( |2.10| , pTT] ). In the case of /'•^■* we simply 



contract eq.( p.lO|) with the polarization vector, S; to calculate we contract eq. 



with a light-like vector, n, which we choose orthogonal to the polarization vector in order 
to put the trace terms in eq.(p.7|) to zero. 



= 0, n-S = 0. (3.37) 



14 



Passing to Euclidean fields and vector components we have 

2M341k = -Sp{pS\i,^Jp^i^\pS), (3.38) 
4M34uark = n^npS^{pS\i,^o.Fp^ij\pS). (3.39) 

To evaluate this matrix element we need the Fourier transform of the dual of the field 
strength of one / and J, eq.( |3.12| ), which is given by 

n2 



= ±p'Gik) ( ^5,. + ^5,, - U,,6,)i , (3.40) 



g{k) = 327r 



i + i)A-„W + (M)A'.W 



t4 



t = kp, (3.41) 



where Kn{t) are modified Bessel functions of the second kind. Inserting this in eq. (|3.34|) 
adding the contributions from J's and /'s, we obtain 

d^k g{k)F{p)F{p - k) 



4M341k = -2M'h{p), 



and 

A /\/r3, 

quark 

. {2-KY{p-kY + M^F\p-k)'i\ k'^p'^ ^ 

where one should set p^ = — after evaluating the integrals. One observes that /quirk 
and rfquark different order in the parameter Mp. For small p both integrals Ii{p) 

and hip) are proportional to p'^p'^, and thus 

41k ~ (Mpf, 

<lk ~ {M pY log Mp (3.44) 

for Mp 0. Numerically one finds, evaluating the integrals for small p, 

2 

hip) = -- P^P^, 

4uLk = ^ = 0.6 (3.45) 



in the limit Mp 0, and 



I^{p) = -0.066 pV, 

y{2) 
^quark 
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^Sark = -0.011 (3.46) 



at the phenomenological value Mp = 0.58. The different parametric order of the twist-4 
and twist-3 matrix elements reflects itself in the numerical values: c^^uark more than an 

order of magnitude smaller than /quark. 

Similarly, we can apply the general formula eq. (|3.34|) to the case of the spin-inde- 
pendent twist-4 operator contributing to power corrections to the Gross-Llewellyn-Smith 
sum rule, eq.( p.l3|) . Contracting eq.( p.l3| ) with the 4-momentum vector and passing to 



Euclidean fields and vector components we have 

2M^41k = P/3|E(pl^7a75^/3a^b). (3.47) 



spm 



Evaluating the matrix element with the help of eq. ( ^.34| ) and summing over spins we 
obtain 

SM^c^lk = -3M2/i(p). (3.48) 

Thus, 



c. 



(2) 

quark 



(Mp)°, (3.49) 



just as the spin-dependent twist-4 matrix element, /quark; eq-( P-44| ), and numerically in 
the hmit Mp 

41k = -I = -0.6. (3.50) 

Finally, we note that an analogous calculation of the forward matrix elements of the 
operators 

{p\^p-faFf3a^\p) , {pS\ij-fa75Ff3ai^\pS), (3.51) 

gives identically zero, as it should be on grounds of the QCD equations of motion. From 
this we conclude that the method of effective operators preserves a principal feature of 
QCD: matrix elements of operators which are zero in QCD due to the QCD equations 
of motion are automatically zero in the effective theory. This remarkable property of 
effective operators will be investigated in more detail in Section p.4| . 

The results obtained here for the quark matrix elements of higher-twist operators 
can be rephrased in terms of local operators, whose matrix elements reproduce those 
of the non-local effective 4-fermionic operators at quark level. At momenta of order 
M <^ p"^ this "localization" is justified by the fact that the effects of the non-locality of 
the instanton-induced vertices contribute only in higher orders of Mp. Since the quark 
loop integrals in eqs. (|3.42| , p.43| , p.48| ) are proportional to the external momentum p^. 



the local operator which should be valid for general quark momenta p of order M, not 
necessarily on the mass shell p^ = — M^, must contain derivatives. The twist-4 operator, 
eq.( |2.10| ), we can thus represent as (in the Euclidean theory) 

i^laFpai^ /quLk ^7/3755V- (3-52) 
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Up to the derivative the equivalent local operator is the axial current operator, multi- 
plied by the coefficient /quark- For the twist-3 operator, eq.( p.7|) , the corresponding local 
operator is 

^ (7^^/37 + 7/3^^7) -dqLk'^{9ad[i3l^]+df3dia-f^])l5ip- (3.53) 



For the spin-independent twist-4 operator, eq. ( P.13| ), finally, the equivalent local operator 
is related to the vector current, 

V'7a75^/3a^ ^ cfu^^rk ^7/35^- (3-54) 

When going to nucleon, for various reasons, a strictly local representation of the higher- 
twist operators is no longer possible. In the nucleon one finds "sea" quarks with momenta 
of order of the cutoff, p~^, so one must take into account the full momentum dependence 
of the instanton-induced vertices. Furthermore, there are interaction contributions to the 
matrix element in which more than one quark connects to the many-fermionic operator. 
The local operators are nevertheless useful as a replacement for the non-local effective 
operator for "valence" quarks whose momenta are restricted to values ^ p^^ by the 
bound-state wave function of the nucleon. Also, the parametric order in Mp which is 
encoded in the constants /quark; '^quLk ^^cl Cq^ark "^ill carry over to the nucleon matrix 
elements. In Section |4.4] we shall see to which extent the relation between /*^^^ and the 
axial vector coupling constant, qa, as well as between c*^^^ and the vector coupling constant, 
which holds for on-shell quark matrix elements, is valid also for the nucleon. 

3.4 The QCD equation of motion 

The set of local operators which are used in the operator product expansion at higher-twist 
level is overcomplete. Some higher-twist operators can be reduced to other operators, 
others are identically zero, by the QCD equations of motion. For instance, the twist-4 
operator 

ij{x)yij{x) = V^(x)7a (^9, - 2yA^(x)^ V^(x) (3.55) 

should have zero matrix element in physical states, 

{P\iP^fi)\P) = 0. (3.56) 
Furthermore, using the fact that in QCD 



y 



F^ai^) = ^[V/3,VJ, (3.57) 



one may show that 

{P\^P^^Ff,^^P\P) = 0. (3.58) 

When passing from QCD to the effective theory it is important that these general proper- 
ties are preserved by the corresponding effective operators. We already noted in Section 
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373| that the matrix element eg. ( |3.58| ) in quark states is identically zero. We now want to 
show in a more general context that eq. (|3.56|) holds in the effective theory, when the gauge 
field in the QCD operator, eq.( |3.55| ), is replaced by the corresponding effective operator. 

Consider the gauge field part of the QCD operator eq.( p.55|) . Since we are interested 
in the forward matrix element we may integrate the operator over the 4-volume, 



d X 



- A" 

'4'{x)la^^{x)Al{x) 



(3.59) 



The effective operator corresponding to this operator is given by eq.( p.21|) p| 



d X 



.(3.60) 



Here the vertices {YA±)a are given by the general formula eq.( |3.18| ), in which one substi- 
tutes the gauge potential of one / (/), eq. (|3.13 ), the Fourier transform of which is given 
by 



/ 



d xe' 



-ik-x 



2l_ 

(x^ + p2)x' 



:X,,6 



(3.61) 



n{k) 



in 



2 r 



t 



kp. 



(3.62) 



Consider a correlation function (connected average) of the effective operator eq. ( 3.601 ) 
with external fermion fields. The connected average consists of two parts, in which the 
external fermion fields connect either to the 2iVj— fermionic vertices, (Ya^)'^, or to the 
factors 7l±: 



Hy') ''Oa" i^Hy) 



il){y') I d^x 
d'^x 



-t^P\xh,^^ix)iYA^)lix) 



^\y)) 



I c 



(7^+) 



A" 



i^^Kx)^,-^lj{x){YA+r,{x) 



N+U{y')'R+^\y) 



+ 



(3.63) 



By explicit calculation, using eqs.( p.2(^ , p.62|) and the fact that the form factor F{k) 
is proportional to the Fourier transform of the instanton zero mode, which implies the 
relation 



d^k , . in — k) ■ k 2z ^, , 



(3.64) 



^Since the instanton-induced vertices {Ya±)'^{x) can have a non-zero vacuum expectation value we 
must use here the full form of the effective operator, including the factors TZ±. 
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one shows that up to terms of order {MpY the insertion of the (2A^/ + 2)-fermionic vertex 
into the correlation function of fermion fields is equivalent to an insertion of a 't Hooft 
vertex, eq.( 



2U{y')Y^ij\y). 



(3.65) 



The factor of 2 appears because, due to the color structure, there are 2Nf different ways 
to contract the {2Nf + 2)-fermionic operator on the right-hand side with the external 
quark fields, but only Nf for the 2A^j— fermionic vertices Y±. On the other hand, the 
vacuum expectation value of the (2iVj + 2)-fermionic vertex appearing in the second term 
of eq.( p.63| ) is 



#t(a;)^^^^(x) {YaMx) 



(3.66) 



by virtue of the self-consistency condition defining the saddle point. The connected 



average of the factors TZ± with external fermion fields is given by eq. (|3.24|) , which is a 
consequence of their definition as the "inverse" of the 't Hooft vertices. Putting these 
results together, one observes that the sum of the two terms of eq.(|3.63D combines to give 



^{y') 'W ^\y) 



^{y')Y±^\y) 



(3.67) 



We have thus shown that, to order {MpY, the gluonic part of the QCD equations of 
motion, eq.( 3.56| ), corresponds to an insertion of the 't Hooft vertex in the effective theory. 
The operator corresponding to the QCD equations of motion thus vanishes by virtue of 
the equations of motion of the effective theory, a very gratifying result. This serves as 
additional confirmation that the approximations made in deriving the effective action and 
the effective operators are mutually consistent. 



3.5 Non-leading vs. leading twist 



In section ^]3| it was shown that matrix elements of gluonic operators of twist 4 are of 
order unity in the packing fraction of the instanton medium, while those of twist 3 are 
suppressed. We now want to investigate the role of instantons in operators of twist 2, 
which determine the scaling (non-power suppressed) part of nucleon structure functions. 
In particular, using the method of effective operators, we want to see why — and under 
which conditions — it is justified to drop the gluon field in twist-2 QCD operators when 
passing to the effective chiral theory. This approximation of "quarks-antiquarks only" is 
used in all calculations of leading-twist structure functions of the nucleon in the effective 
chiral theory |16|. 



Consider the twist-2 QCD operator of spin n + 1, eg. ( p. 41) , which determines the ra'th 
moment of the polarized structure function, gi (for simplicity let us consider one flavor, 
Nf = 1). With the explicit form of the covariant derivative, 

= ^,-^^A;,ix), (3.68) 
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the operator can be written as a part containing only explicit derivatives, plus terms 
involving powers of the gauge field|], 

^7{a75V;3, . . . V^„}V^ = V^7{a75<9/3i • • • d 

n 

i ^ 
+ terms A^ . . (3.69) 

When passing to the effective chiral theory, the pure derivative part becomes simply the 
corresponding operator in terms of fields of the effective theory, while the terms involving 
the gauge field have to be replaced by the corresponding effective operators; they are of 
the form of the general operator eg. ( ^.151 ). Let us consider the matrix elements of eq.( |3.69| ) 



in quark states of the effective theory. The matrix element of the pure derivative part is 
trivial and contributes precisely unity to the moment a^'^\ On the other hand, using the 
general formulas of sections |3]^ and |3]^ one easily shows that the matrix element of the 
terms involving the gauge field is of order (Mp)^. Thus, 



(n) 



quark 



1 + O [{Mpf) . (3.70) 



The contribution of the gauge field is parametrically suppressed relative to the pure deriva- 
tives. In terms of the structure function this implies 

^l(x)quark = \b{x - \) ^ O [{Mpf) . (3.71) 

To leading order in Mp the quark of the effective chiral theory is "elementary" , it structure 
in the instanton vacuum emerges only at level (Mp)^. 

At the same time, one finds that the twist-2 QCD operators determining the moments 
of the polarized gluon distribution, 

{PS\F{^pDp, . . . Dp^_,Fl^\PS) - traces 
= 2M;v&^"^ 5{aP/3i...P/3„} -traces, (3.72) 

are identically zero when one substitutes the field of one instanton (the same applies to the 
unpolarized gluon distribution.) The reason is the 0(4) symmetry of the instanton field: 
after integration over instanton coordinates one can build the (3i . . . (3n tensor only out 
of Kronecker symbols, so it is impossible to get it traceless. Thus the effective operator 
for eg. ( |3.72|) is zero. However, the operator in eq.( |3.72|) becomes non-zero when taking 
into account two-instanton contributions, which are parametrically suppressed by a factor 
{p/RY oc [Mpf. We conclude that 

= 0({Mpf). (3.73) 

To summarize, we see that in twist-2 operators the effect of the instanton field is 
always of order {MpY. When working to leading order in Mp, it is consistent to keep 

''We remind that we have in mind here the gauge in which the instanton field takes the form eq. ( |3.13| ) , 
the so-called singular gauge. 
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only the pure derivative part in the twist-2 quark operators, eq.( p.69| ), and to replace the 



operators for the twist-2 gluon distribution by zero. This is the approximation employed 
in the calculation of the twist-2 part of the nucleon structure functions in refs. []T3|, |TB|, [T^. 

However, as we have seen above, instantons can make contributions of order unity 
in higher-twist operators. In the cases considered in section the first non-suppressed 



instanton contributions appeared at twist 4. There can of course also be 0((Mp)^) con- 
tributions to operators of twist 6 and higher. 

The difference between the instanton contribution to operators of leading and non- 
leading twist can be seen even more clearly in the case of unpolarized structure functions. 
The second moment of the structure function F2, M^'^\ is given by the matrix element of 
the twist-2 spin-2 operator, 

(P|V'7{mV^2}V^I^) -traces = M(2)p^^^ P^^^ - traces. (3.74) 

In analogy to eq. (|3.70|) , at quark level the pure derivative term makes unity contribution 
to M^'^\ while the contribution of the gauge field is suppressed by {MpY". 

= 1 + {{Mpf) (3.75) 

This is consistent with the energy-momentum sum rule, since the momentum carried by 
gluons is 0((Mp)^), cf. eq. (|3.73|) . On the other hand, if one does not project on twist 2 



but contracts the indices in eq. (|3.74| ) one obtains the twist-4 matrix element 



(P|^y^|P) = CM^, (3.76) 



where C = by virtue of the equations of motion. As was is shown in section ^.4| , this zero 
comes because an order unity contribution to C from the instanton field exactly cancels 
the pure derivative term in the Dirac operator. 



4 Nucleon matrix elements of higher— twist operators 
4.1 The nucleon in the effective chiral theory 

An immediate application of the effective chiral theory derived from the instanton vacuum 



is the chiral quark soliton model of the nucleon |13|. In this section we briefly recall the 
main characteristics of this description. 

In the effective chiral theory the nucleon is described in the large- A^^c limit by a classical 
pion field, Uc{'x)- The quarks single-particle wave functions are determined from the Dirac 
equation in the external pion field, 

iJ$„(x) = P„$„(x), (4.1) 



H = 74Mfc-Mt/^«(x)], (4.2) 
W'^ = exp[m'^(x)r"75] = ^^^^-^U+{x) +^^—-^U\x). (4.3) 
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The spectrum of the one-particle Hamiltonian, H, contains the upper and lower Dirac 
continua (distorted by the presence of the external pion field), as well as a discrete bound- 
state level. We denote the energy of this discrete level by Ei^^, —M < Ei^^ < M. This 
level must be occupied by Nc quarks to have a state of unit baryon number. The saddle 
point pion field is determined by minimizing the static energy, which includes the energy 
of the discrete level as well as the aggregate energy of the negative Dirac continuum. It 
is of "hedgehog" form. 



f/c(x) = exp [i7r"(x)r''] = exp 



F(r] 



(4.4) 



where P(r) is called the profile function. A variational approximation to the soliton profile 
in the form 

P(r) = 2 arctan (4.5) 

with ro ~ 1.0/M gives a very reasonable description of a variety of static nucleon observ- 
ables [|I3]. 

The minimum of the energy is degenerate with respect to translations of the soliton in 
space and to rotations of the soliton field in ordinary and isospin space. For the hedgehog 
field eq. (|4.4]) the two rotations are equivalent. Quantizing slow rotations of the saddle- 
point pion field gives rise to the quantum numbers of the nucleon: its spin and isospin 



components ||28|, [T^. In order to take into account the translational and rotational zero 
modes one has to perform a rotation of the soliton field in flavor space and to shift its 
center, 

t/(x) = RUc{^-X)R\ (4.6) 

the same for the quark eigenfunctions, 

$„(x) ^ i?<l>„(x-X), (4.7) 

and to make a projection to a concrete nucleon state under consideration. The projection 
on a nucleon state with given momentum P is obtained by integrating over all shifts X 
of the soliton, 

(P'l . . . I P) = I e^(P'-P)-x .... (4.8) 



The projection on a nucleon with given spin (S'3) and isospin (T3) components is obtained 
by integrating over all spin-isospin rotations R of the soliton, 

= T, S3, Tal . . . |S = T, ^3, T3) = jdR 4>\!n^R) ■ ■ ■ <Pl7T,iR), (4.9) 

where (/)|^(i?) is the rotational wave function of the nucleon given by the Wigner finite- 
rotation matrix ||13|| : 

4;^^{R) = v/25TT(-l)™I)!g53(^)- (4-10) 

The four nucleon states have S = T = 1/2, with 5*3, T3 = ±1/2. It is implied that dR 
in eq. (|4.9|) is the Haar measure normalized to unity. (In the following we shall omit the 
superscript S = T.) 
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4.2 Nucleon matrix elements of effective gluon operators 



We now turn to the calculation of nucleon matrix elements of the effective operators 
which arise in the description of higher-twist corrections. In this subsection we derive 
the expressions for the nucleon matrix element of the general gluonic operator, eq. (|3.15 ), 



which was considered in Sections 3.2 and 3.3 



Oai...ar(3i...Psi^) 



where we have introduced now a flavor matrix, t, with i = for the fiavor-nonsinglet 
and t = 1 for the singlet case. 

To compute the nucleon matrix element we use the bosonized form of the effective 
operator, eq.( p.32| ), where the pion field is to be replaced by the soliton field, eq.( [4.4|) : 



"0"oi...Q,,,/3i.../3^(x) 



X 



A" 



'+,fivPx...l3s 



fJ.uf3l...Ps 5 



1 + 75 1-75 



f/4 



(4.11) 



(For simplicity we have neglected for the moment the nonlocality of the instanton-induced 
vertices, i.e., we have set the form factor F{k) = 1. The form factors will be inserted again 
later.) This operator has the form of a product of two color-octet currents, similar to 
one-gluon exchange between quarks, with the functions J^±,^vp^...p^ playing the role of the 
gluon propagator. The calculation of the nucleon matrix element of the effective operator 
is thus analogous to that of gluon exchange corrections to the nucleon mass One 
expands the quark fields in eq. (|3.32|) in the basis of single-particle wave functions, eq.( [4.1|) , 
and computes the matrix element, summing over all occupied quark single-particle levels 
including the bound-state level and the negative Dirac continuum. In addition, one has to 
perform a rotation and shift of the center of the pion field and the quark wave functions, 
eqs. (^I6| , ^]3), and integrate over collective wave functions to project on a nucleon state 
with given spin-isospin quantum numbers. One obtains 



(P — 0, T3, S'sl "0"q,^...q,^/3^.../3JP — 0,T3,S'3) 
= 2M^ / dRct>\^T,{R)4>s,TAR) -M«i...«./3i.../3.(i?), 



(4.12) 



M 



a^...arf3i.../3s{R) 



I d'^x I d'^z J^+,f,uf3i...f3s{x - z) Xl'exp [(E„ - Em){x^ - z^\ 



X 

+ (+ 



<l>m(x) <l>J„(z)74 a 



1+75 



MUM $Jzl 



(4.13) 
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where 



(4.14) 



m,,n 



n m 
occ. non— occ. 



n m 
non— occ. occ. 



Here we have used that for the forward matrix element, P = P' = 0, the integral over 
translations is trivial. Due to translational invariance it is equivalent to integrating the 



coordinate of the operator, x, eq.( |3.32| ), over the 3-volume. To obtain a more symmetric 
form, we have integrated over the Euclidean 4-volume and divided by the Euclidean time 
interval, T. In the following we shall drop the labels P = 0. 

A factor of resulting from the sum over color indices combines with the factor 

l/Nc in the effective operator eq. (|3.32|) , to give an overall factor of Nc/2. 

To proceed with the evaluation of eq.( [4.13| ) one has to perform the integral over soliton 
rotations, considering separately flavor singlet and nonsinglet operators. We note that 
the factor $J„ . . . U± . . . in eq.( [4.13| ) does not depend on the rotation matrices, as it 



should be. This factor represents the "gluonic" part of the operator eq.( 4.11 ), and is thus 
by definition a flavor singlet. 

One easily sees that for spin-dependent matrix elements such as S'^^ and /'■^\ eqs. (|2.7| , 
|2.10| ), the flavor-nonsinglet part appears in the leading order of the 1/iVc-expansion. In 
this case in eq.( [4.13| ) t = r^, and the integral over rotations is performed using 



■tr 



R^t'^Rt^ 



(4.15) 



dR(Pl^j,.^{R)Dsa{R)(pSsTs{R) 

The non-singlet matrix element thus becomes 

/rri Q I U/^JlNS \rri Q \ 

\J35'J3| <^ Qi...ar/3i.../3j-'35 '-'3/spm-dep. = 



-2M^(2T3)^A^«^...„^^^...^^, 



(4.16) 



•'^'ai...Qr/3i.../3, 



X 

+ (+ 



$^(x) $J„(z)74a, 



1+75 



MU+{-L)^n{^) 



(4.17) 



Here Sa (a = 1,2,3) are the components of the nucleon polarization vector in the rest 
frame, 



S 



(0,S), 



(0,0,1). 



(4.18) 



For spin-independent matrix elements such as c*^^\ eqs. (|2.13| ), on the other hand, the 
flavor singlet part is leading in N^- In the flavor singlet case t = 1, the orientation matrices 
i?, i?^ are contracted, and the integral over rotations becomes trivial. 



dRct>i^^.^{R)<Ps,T,{R) 



1. 



(4.19) 
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One obtains 



Ma,...a,.f3^.../3. = (4.17) with ^ 1. (4.20) 

The isosinglet spin-dependent and isovector spin-independent matrix elements are 
zero in the leading order of the 1/A^c^expansion. They appear only after including ro- 
tational corrections to eq.( |4.13| ), i.e., integrating over time-dependent rotations and ex- 
panding in the angular frequency, which is of order l/N^. 

The A'^c-ordering of the matrix elements of twist-3 and 4 operators is analogous to 
the one of the twist-2 quark distribution functions, where the isovector-polarized and 
isosinglet-unpolarized distributions appear in the leading order of the l/A^^c-expansion. 
The connection between spin and isospin is a general consequence of the hedgehog sym- 
metry of the pion field, eq. ( [4.4D . For reference, we have listed order of the various matrix 
elements in in table |1|. 

We have obtained a closed expression for the matrix element of the effective operator 
corresponding to a general gluonic operator eq.( p.l5| ), in the form of a sum over quark 



single-particle levels in the pion field of the soliton. This representation is useful for 
studying the A'^c-order of different spin-isospin combinations. Moreover, eq. ( |4.13| ) could 
serve as a starting point for numerical calculations of the matrix elements, using the 
quark wave functions determined by numerical diagonalization of the Dirac Hamiltonian, 



eq. ([4.1| ), in a finite box |T^, For further theoretical analysis, however, it is convenient 
to pass to the more familiar language of Feynman diagrams. Substituting in eq.( [4.13| ) the 
spectral representation of the quark Green function in the background pion field, 

G(x, y) = e{x4 - 2/4) J2 6xp [-Enix4 - 1/4)] <l>„(x)$J^(y)ii74 

E„>0 

- ^(i/4-X4) 51 exp[-E„(x4-?/4)]$n(x)$],(y)i74, (4.21) 

E„<0 



and rearranging the terms in the sum over quark levels, one can show that eq. (|4.13| ) 
is equal to the sum of three terms which correspond to the Feynman diagrams for the 
3-point correlation function shown in Fig.p|, 



(^3, S^l IT3, 5'3)spin-dep. 

-2M^(2T3)% \MI,^ + Ml,,) + M1,,,)\ ^ „ . (4.22) 



3 



ai...Or/3i---/9s 



We write them down for the spin-dependent isovector matrix element, eq.( ^4.17| ); the 



corresponding expressions for the spin-independent isoscalar matrix elements are obtained 
by the trivial substitution eq.( [4.20|) : 



(i) A "valence quark" contribution in which the effective operator connects to one 
valence quark line: 

-'^'Qi...a,/3i.../3,,(i) 
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X I exp [Eiev (X4 - 2:4 ) ] ^\^^ (x) 
+ exp [-Eiev(x4 - Z^)] $L(Z 



74rai...a,r°G(x,^;)74a 



1+75 



MUjz] 



$lcv(z) 



1 + 75 



+ +--. 



$lev(x) 

(4.23) 



This contribution has the form of a quark self-energy diagram evaluated between bound- 
state wave functions. This contribution is similar to the Lamb shift in atomic physics. 

(ii) A "sea quark" contribution, in which the effective operator is not connected to the 
valence quark lines but to a closed quark loop in the background pion field. In momentum 
representation it can be written as 



■^ai...arfSi...l3,,{ii) 



Nr 



d^h 



(277)4 J (27r)4 J (27r)4 J {2n)^ J {2tt) 



d% 



d% 



X (27r)^5(^) {q-h + h + h- k2)2Ti6{q^) 
X ^+,^./3,.../3,(A;2 - A;i)F(/i)F(/2)tr 
+ (+--), 



^(^2, h)T-T^„„^Mki. /i)a^.i^Mf/+(q) 



(4.24) 



where the quark Green function in the background pion field in momentum representation 
is defined as 



G{x,z) 



d^k 



dH 



{2n)^J (27r) 



^G{k,l), 



and 



d^xe-''^ '' [f/±(x) - 1] 



(4.25) 



(4.26) 



is the Fourier transform of the hedgehog pion field, eq.( [4.4|) . In eq.( [4.24|) we have rein- 
stated the form factors F{k) which cut the quark loop momentum. 

(Hi) An "interaction" contribution, in which the effective operator connects two dif- 
ferent valence quark lines: 



■^ai...ar(3i...l3s,{m) 



2T 



d X d zT. 



+,fiuf3i...p. 



{x — z) 



X 



*L(x)74T"r„,...„,$lev(x) 



$L (Z)74a^. Mf/+ (Z) $iev (Z) 



+ (+--)• 



(4.27) 



This contribution describes interactions of pairs of valence quarks in the nucleon mediated 
by the gluonic operator. It is similar to the hyperfine splitting of baryon masses due to 
one-gluon exchange 
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We emphasize that the expressions derived in this section apply to the leading order 
of the 1/A^c^expansion. For matrix elements which are zero in leading order of the l/N^- 
expansion one has to include rotational corrections. 



4.3 Higher— twist matrix elements in the naive quark model 

Before evaluating the nucleon matrix elements of higher-twist operators in the chiral 
soliton model it is useful to consider a much simpler picture of the nucleon, a naive 
constituent quark model. Assuming N^. free constituent quarks in a spin-isospin state 
described by the non-relativistic quark model, one can easily construct the nucleon matrix 
elements from the quark matrix elements derived in section |3.3| . We wish to caution the 
reader that, as will be seen in the following section, the effects of the binding of the quarks 
are essential, and the quark model results turn out to be unreahstic; nevertheless, it is 
useful to have at hands these simple estimates. 

In section it was seen that the quark matrix element /quai.^ '^^'^ ^e reproduced by 
the local operator eq.( p.52| ). For on-shell quarks the derivative in this operator simply 
becomes (9^ = M^, hence the operator reduces to the axial current. This immediately 
tells us that in the naive quark model f^"^^ is proportional to the nucleon axial coupling 
constant (for both flavor-nonsinglet and singlet) 

f(2) _ A2) (3) .(2) _ M2 (2) (0) 

Ins — /quark 9a : JS — /quark 9 A ■ l^-^oj 

(3) 

The nucleon isovector axial coupling constant in the large- A^c quark model is g\ = 
{Nc + 2)/3 = 0{Nc), while the isosinglet coupling constant g^^^ is 0{N^), so the relative 

(2) (2} 

order in 1/Nc of f^g and is in agreement with the general results of the last section. 
With /quLk given by eq.( p.45| ), the phenomenological value g^^^ = 1.25, and a value of 



g\ = 0.3 obtained from an analysis of polarized nucleon structure functions of [^, this 
would mean positive values for f^^l = 0.75MVM^, = O.ISMVM^. 

Similarly, for the spin-dependent twist-3 matrix element d^'^\ eq. (p^ ), one obtains in 
the naive quark model 

(2) _ (2) (3) (2) _ (2) (0) 

~ j^2 '^quark 9a y "S ~ '^quark 9a 5 1^4. Zyj 

which with eq.dCTj) would imply d^^l = -0.013 M^/Mjj and = -3-3 x lO^^MVMlf. 

The spin-independent matrix element c^'^\ eq.( p.l3| ), is at quark level related to the 
matrix element of the vector current-type operator, eq.( p.54| ). Thus, in analogy to the 
spin-dependent case one finds that in the naive quark model Cg is related to the nucleon 

f 2") 

vector charge (baryon number), which is Nc^ while is proportional to the nucleon 
isospin, which is of order 0{N^), 

(2) _ (2) ^ (2) _ (2) .^3Q^ 

~ Suark ^^c, <^NS ~ ^uark' {fi-ClU) 

The quark model estimates are useful since they illustrate the Ai'c-ordering of the 
different matrix elements and their parametric magnitude in Mp. We shall see, however. 
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that the effects of the binding of the quarks in the nucleon cannot be neglected because 
of the momentum dependence (non-locahty) of the effective operators. 



4.4 The leading higher— twist matrix elements in the large— A^^ 
limit 

We now evaluate the leading twist-3 and 4 matrix elements in the large-A^^c limit in 
the chiral soliton model of the nucleon. As was shown in section the leading matrix 
elements are the flavor-nonsinglet spin-dependent matrix elements rf^g and fj^l, eqs.( |2.7] , 
|2.10|) , and the flavor-singlet spin-independent matrix element Cg \ eq.( |2.13|) . 

In addition to the 1/iVc-classification we must keep in mind that our theory is based 
on the smallness of the parameter Mp. In section we have seen that at quark level 
the twist-4 matrix elements f^"^^ and c^^-* are parametrically of order unity, while S"^^ is 
suppressed by a factor of {MpY . We expect the parametric order of the nucleon matrix 
elements to be the same as that of the quark matrix elements; this is indeed the case, as 
will become clear below. We shall therefore now consider separately the parametrically 
large (twist-4) and small (twist-3) matrix elements. For quick reference we have listed 
the parametric order of the various matrix elements in table 1^. 

Let us first consider the spin-dependent twist-4 matrix element f}^g, eq.( |2.1CI|) . A 



parametrically large contribution can come only from "divergent" diagrams in Fig.g, i.e. 
diagrams with characteristic momenta of the order of the cutoff, p~^. The diagrams con- 
taining divergences are the valence quark contribution (z) and the sea quark contribution 
(ii), while in the interaction diagrams {Hi) all momenta are cut by the bound-state level 
wave function. We shall therefore concentrate on extracting the leading contribution in 
Mp from diagrams {i) and (u). 

The valence quark contribution diagram (i), eq. ( [4.23| ), contains a quark "self-energy" 



diagram of the same type as appears in the on-shell quark matrix element, see Fig.Q. 
In the nucleon case, of course, the momenta of the external quark lines are off-shell, as 
determined by the bound-state level wave function. Since these momenta are of the order 
of the quark mass, M, it is sufficient to compute the self-energy subdiagram for momenta 
of order M — one does not need to know its full momentum dependence up to momenta 
of order p"^ . In other words, one may compute the contribution {%) to the nucleon matrix 
element f^^\ using the local operators eq.( p.52| ); corrections to this local approximation 
are of higher order in Mp. We thus estimate the contribution eq. ( [4.23| ) to f'^\ as 

.(2) _ ^ .(2) ~(3) 

JNS,{i) ~ f./r2 i quark i/A,lev? V^-"^^) 



N 



where g^^Xev denotes the matrix element of the isovector part of the operator eq.( |3.52| ) 
between the bound-state wave function^, 

~9%. = ^/ci'a;<l>L(x)74[r"7a75(i?L + ^')]<^'icv(x), (4.32) 



*The explicit form of the bound-state wave function can be found in ref.jn 
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(3) 

which is analogous to the level contribution to the axial coupling constant, (/Viev) but con- 



tains the derivative operator. The effect of the latter are dramatic. Evaluating eq. (|4.32| ) 



with the variational soliton profile eq.( [4.5|) with tq = 1.0/M we find g^lev — —0.75, while 



the level contribution to the axial coupling constant is g\[^^ = 0.70. Note that for a 



free quark, as was assumed in the naive quark model estimate of section [4.3|, one would 



have ^^lev = (^Aiev Thus, the binding of the quarks by the pion field reverses the sign 
of the matrix element of the derivative operator eq. (|3.52| ). This fact is not dependent on 
the precise value of the parameter ro which characterizes the radius of the classical pion 
field and thus the binding energy of the bound-state level; one finds that gA,\ev becomes 
positive only for unphysical values tq < 0.5/M where the level is practically unbound. At 
the physical value tq = 1.0/M we thus obtain, cf. eq.( |3.45|) . 



f(2) 

JNS, (i) 



-0.45 X 



(4.33) 



In the sea quark contribution, diagram (ii), eq. (|4.24| ), on the other hand, all quark 
momenta are of the order of the cutoff, p~^, hence we must take into account the full 
momentum dependence (non-locality) of the vertices in the effective operator. To render 
an evaluation of this diagram feasible one can make use of an expansion of the quark 
propagators in the background pion field whose characteristic momentum is of order M. 
The expansion can be performed in such a way that it becomes exact in three limiting 
cases: i) small momenta, ii) large momenta, and Hi) any momenta and small pion fields. 
This is the so-called interpolation formula of ref . [l^ . Applying this expansion to eq. (|4.24| ) 
we obtain 



JNS, (ii) 



S'q 
(2^ 



/(q) qa tr 



-q)t>+(q) 



(4.34) 



where J(q) is given by a two-loop integrals with quark propagators in the trivial back- 
ground field (?7±(x) = 1), 



X 



G{k)F{k2)F{k^)F\k^ + q) 



(27r)4 J (27r)4 {kl + M2F4) {kl + M'^F^) [{ki + g)2 + W^F^)] 
'^N{h,k2,q)+^-N{h,k2,Q) 



(4.35) 



N{ki,k2,q) = - i-ki ■ {ki + q) + ^-^ ^ \ + k-{ki + q), 

k = k2 — ki. 

Evaluating the integrals numerically at Mp = 0.58 we obtain 

M2 



A2) 
JNS, (ii) 



-0.27 X 



(4.36) 



(4.37) 



Adding the contributions (i) and (ii), eq. ([4.33|) and eq. (|4.37| ), using the standard value 
for the quark mass, M = 350 MeV, and Mn = 940 MeV, we obtain our estimate for the 
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nucleoli matrix element, 



fSl = -0-72 = -0.10. (4.38) 



The interaction contribution, diagrams (iii), eq.( |4.27| ), is parametrically of order 
{MpY. Explicit caclulation of the integral in eq. ([4.27|) shows that this contribution to 
f)y[g is numerically an order of magnitude smaller than the one from the "divergent" dia- 
grams {%) and (m) and can safely be neglected. (We note, however, that in parametrically 
suppressed matrix elements such as S-"^^ the interaction contribution [iii] is of the same 
order as (z) and {ii) and must be included.) 

It is interesting to consider the theoretical limit of large soliton size, the "skyrmion 
limit", which may be regarded as the opposite of the naive quark model limit|. In this 
limit the valence quarks are strongly bound, i.e., the energy of the discrete bound-state 
level approaches the lower continuum. In this case diagrams [ii) account for the entire 
contribution. Furthermore, one can neglect the momentum dependence of the form fac- 
tor I{q) in eq. ( |4.34| ), in which case the function of the pion field in eq.( [4.34| ) becomes 



proportional to the axial coupling constant, which for large solitons is given by 



--f / d^x tr l-iT^U^daUi^)] , (4.39) 



where is the pion decay constant, which is defined by a logarithmically divergent 
integral, 

Fl = 4iV,M2/-^ ^ (4.40) 

(Numerically, = 93MeV for Mp = 0.58.) Thus, in the limit of large soliton size we 
obtain the relation 

^ -Sf^).J' ^ -0.53 x^«- ,4.41) 



similar to the naive quark model relation eq.( |4.28| ), but with a negative coefficient, since 
/(O) is positive, /(O) = 0.95 x 10~^. Quantitatively, with g^^^ = 1.25 this relation works 
well when compared to the result for physical soliton sizes, eq.( [4.38| ). 



The calculation of the spin-independent twist-4 matrix element c^g\ eq. (|2.13|) , pro- 
ceeds in much the same way as that of f^g. This matrix element is also of order unity in 
Mp, so the dominant contributions come from diagrams (i) and (ii). The valence quark 
contribution can again be estimated by "localizing" the self-energy subdiagram in (i): 

(2) _ (2) (0) .,,2) 



^For hadronic observables the two limits were investigated in ref.f 
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where gyl^^ is the isosinglet matrix element of the local operator eq. ( |3.54| ) between bound- 
state wave functions, 



~{0) 
9vAev 



(4.43) 



Without the derivatives this would be just the isosinglet vector charge, i.e., the baryon 
number, Nc- Numerically we find Qylev = —1-43 x Nc- The sea quark contribution (ii) to 

(2) 

Cg emerges only in the third order of the expansion in pion field momenta (interpolation 
formula), analogous to the expansion of the baryon charge, and is thus strongly suppressed 
compared to the valence quark contribution for physical soliton sizes. Our result for the 
nucleon matrix element c^p is therefore given entirely by eq. ( |4.42| ), 

M2 



.(2) 



2.6 X 



Ml 



0.36. 



(4.44) 



(2) 

Finally, we note that in the limit of large soliton size the nucleon matrix element 
becomes proportional to the winding number of the pion field. 



Q 



247r2 



{U^d,U)iU^d,U)iU^dkU) 



d oc ^ tj r 

which for large solitons coincides with the baryon number [|T^]. One obtains 

M2 



^(2) 



Ml 



where 
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727r^ {MpY ho N,Q, 



g{k) F\k2)F{kr 



(4.45) 



(4.46) 



(27r)4 J (27r)4 [kl + M^F^{k2)r [kf + M^F\ki)] 



N{ki,k2 



N{hM) = + ^ = (4.47) 

One finds that in the limit Mp —>■ 

{Mp)^ho ~ (Mp)°, (4.48) 

in agreement with the fact that Cg^ is parametrically of order (Mp)^. 

Let us now turn to the spin-dependent twist-3 matrix element, d)^s, eq.( |2.7| ). At quark 
level this matrix element is proportional to {MpY. At the present stage of development of 
our theory we cannot make rigorous predictions at level (Mp)^; this would require, among 
other things, to go beyond the one-instanton approximation in the effective operator. 
Nevertheless, it is worthwhile to make a crude estimate of the nucleon matrix element by 
evaluating the one-instanton effective operator in the nucleon. Since for parametrically 
suppressed matrix elements the effects of the non-locality of the effective operator are 
important, we consider the limit of large soliton size where the entire contribution is given 
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by diagrams (ii). To compute c?^^ we contract eq.( |4.24| ) with the "hght-hke" Euchdean 
vector n as in eq. (|3.39|) . This vector must now be taken in the nucleon rest frame, 



n 



U, n 



n- S 



0. 



(4.49) 



(2) 

We then proceed as in the case of fj^g, expanding the quark propagators in derivatives of 
the pion field. One finds that the first non-zero term appears at third order, 



C[U] 



_M2 N, 

niTljSkSa 



{Mpf /i2 C[Ul 



X / d X tr 



(4.50) 

+ ^ 3] 

(4.51) 



Inserting the hedgehog pion field, eq.(O), and averaging over orientations of the polar- 
ization vectors eq. (|4.51| ) becomes equal to 



C[U] 



167r2 



drr"^ 



2 sin P cos P ^/ 2 sin^ P 6 sin'^ P cos P 

P P \ 3 



. (4.52) 



Numerically, C[U] = 0.52 x 10^ for the variational soliton profile of eq. (|4.5| ). In eq.( 4.50| ) 
/i2 is defined as 



'12 



d'^ki f d'^k2 



g{k) F\k2)F{k^f 



(27r)4 J (27r)4 [kl + M^F^{k2)Y [kj + M^F^ik^)] 



rN{ki,k2), (4.53) 



and one finds that 



{Mpfh2 



oc 



Mp 



(4.54) 



for Mp — > 0, in agreement with the fact that at quark level d*^^^ is parametrically sup- 
pressed. Numerically, at Mp = 0.58 one obtains (Mp)^/i2 = 0.43 x 10~^, so that 



(2) 
NS 



-0.1 X 10 



-2 



M2 



(4.55) 



This number should be taken as an order-of-magnitude estimate. 

To summarize, we find that the nucleon matrix elements of higher-twist operators are 
of the same parametric order in Mp as the corresponding matrix elements in quark states. 
(For reference, the parametric order of the different matrix elements is summarized in 
table |I[) Also, their numerical order of magnitude agrees with the quark model estimates. 
However, we see that the binding of the quarks in the nucleon plays a decisive role in the 
twist -4 nucleon matrix elements: the effective operators at quark level are non-local (in 
the simplest form they may be approximated by local operators with derivative (9^), and 
the finite size of the bound state reverses the sign of the matrix element as compared to 
free quarks. We conclude that for a reliable description of twist-4 matrix elements it is 
of crucial importance that the operators and the nucleon state are treated consistently 
within the effective theory. 
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twist 


0{Mp) 


0(iV,) 




4 


1 




l/N, 
1/N^ 




3 


(Mp)2 


S 


1/iVc 
1/N^ 


C(2) 


4 


1 


S 
NS 


1/iVc 



Table 1: A summary of the parametric order of the spin-dependent higher-twist matrix 
elements Z^^) and eqs.(^, ^J^, and the spin-independent matrix elements \ 
eq. (|2.13|) , in 1/Nc and the parameter Mp related to the packing fraction of the instanton 
medium, cf. eq.( \i.T\ ). Details are given in the text. 



5 Results and discussion 



We now want to compare the numerical results for the twist-3 and 4 matrix elements 
with theoretical estimates obtained using other non-perturbative methods, as well as 
with estimates derived from experimental data on polarized nucleon structure functions. 
Numerical values are listed in table 0, where we also quote results of QCD sum rule 

estimates obtained by Ji and 
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calculations by Balitskii et al. and Stein et al. [| 
Unrau using the bag model |]21| , as well as lattice results from Gockeler et al. ||3^ . 

When discussing numerical values we must keep in mind that the theoretical status 
of our numbers presented in table |^ depends on their parametric order, see table |l]. The 



results for quantities of order unity in Mp, 



may be regarded as quantitative 



estimates; the numerical value quoted for the parametrically suppressed twist-3 matrix 
element, d)^si should be taken as a crude estimate. We can say confidently only that 
numerically the twist-3 matrix elements are more than an order of magnitude smaller 
than the twist-4 ones. 

Our estimate for f'^\ argees well with the results of QCD sum rule calculations of 
and [|^, In particular, the large non-singlet value for f^'^\ which was first observed 
in the sum rule calculation of Balitskii et al. [Q, can naturally be explained by the fact 
that the non-singlet part is leading in the l/A'^c^expansion. Also, our result for the spin- 
independent matrix element Cg (corrections to the GLS sum rule) is in good agreement 
with the sum rule calculation of Braun and Kolesnichenko [^. In their calculation the 
fiavor-singlet part was found to be large, which, again, is consistent with the l/N^- 
expansion. 



(2') 

We also note that our result for /^^ differs in sign from the bag model value 



of ref. 21 



As to the twist-3 matrix element, 



(2'] 

^NSi "^oie that our small value 



typical scale of the twist-4 matrix elements, which is M'^/Mjf 



10" 



"small" on the 
are consistent 



with the results of measurements of the polarized structure function g2 ||3^, ^ ^ , where 
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/■(2) 
J NS 






scale/GeV^ 


present work 
Balitskii 33 
Stein 34, 35 
Ji [21] 

Gockeler g 
Braun 24 1 

E142, E143, E154 |37, 38| 
Ji 13 


-0.10 
-0.20 
-0.072 
0.11 

0.10 ±0.28 


~ 10-3 
0.072 
0.072 
0.063 
-0.13 

0.043 ± 0.046 


0.36 
0.37 


~ 0.4 
1 
1 

5 
4 
1 
3 
1 



Table 2: Summary of numerical results for the flavor-nonsinglet spin-dependent twist- 

4 and 3 matrix elements f^g and (ij^g, eqs.(p.7|, p.lO| ), and the flavor-singlet spin- 
as) 



independent twist-4 matrix element , eq.( |2.13| ). Flavor-nonsinglet and singlet parts 



are deflned as in eqs. (|2.3| , |2.4|) . Shown are the results of the present calculation (for details 
see Sections [4.4[), of QCD sum rule calculations by Balitskii et al. 



Stein et al. 



and by Braun and Kolesnichenko p4[, results of a bag model calculation by Ji and Unrau 
as well as lattice calculations by Gockeler et al. 
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Also shown are estimates of 



S^g based on measurements of the structure function g2 by E142, E143, E154 p7 



(for details see the text), and estimates of fl^g from a recent analysis of power corrections 



to gi by Ji and Melnichouk |^ 



S'^^ enters in the non-power suppressed-part, eqs.( p.8| , |2.9|) . (The flavor-nonsinglet value 



we quote in table || has been combined from the E143 value for (P''> and the SLAG 



average for c?^^^ based on measurements of E142, E143 and E154 at a scale of 
3GeV^.) In contrast, the lattice results of seem to overestimate these quantities. 



Finally, we note that using the numerical results of table ^ one can establish the 
Q^-dependence of the Bjorken, Ellis- Jaffe and Gross-Llewellyn-Smith sum rules. For 
this one has to add to eqs.( |2.1| , |2.2|) the logarithmic corrections, which can be found in 
PI]. Let us mention also that our value for f'^\ is within the bounds obtained from a 



phenomenological analysis of power corrections to gi by Ji and Melnichouk ED] 



6 Conclusions and outlook 

In this paper we have reported results of a comprehensive investigation of matrix elements 
of operators of twist 3 and 4 in the instanton vacuum. We have seen that, generally, the 
instanton vacuum provides a clear and consistent picture of the role of non-perturbative 
gluonic degrees of freedom in DIS matrix elements. The crucial element is the small 
parameter Mp oc (p/R)'^ inherent in this picture. In particular, the fact that the QGD 
equations of motion are realized at the level of effective operators provides a check for the 
consistency of this approach. 

A marked distinction between operators of highest twist (twist 4, in our case) and 
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lower twists has emerged. The contribution of the instanton field in operators of twist 2 
and 3 is suppressed by a factor of the instanton packing fraction, (Mp)^. This justifies 
dropping the gluon field in twist-2 operators {i.e., replacing covariant by pure derivatives) 
when computing twist-2 quark and antiquark distributions in the effective chiral theory 
T3| , |T^. We see at present no other approach which provides an equally clear — that is, 



parametric — justification for this "quarks-antiquarks only" approximation which proves 
to be very successful in practice. On the other hand, instantons make an order unity 
contribution in operators of twist 4. The most striking example is the twist-4 operator 
of Section |3.4| which is zero by the QCD equations of motion: in the instanton vacuum it 
vanishes because of an (Mp)° instanton contribution. 

The deeper reason why instantons contribute only to the highest twist considered here 
is the 0(4) symmetry of the instanton, as a result of which the instanton can contribute 
only to the lowest-spin projection of a given tensor operator. To which twist this cor- 
responds depends on the situation considered. Thus, rather than speak about twist one 
should say that instantons prefer lowest spin. 

As a result of the parametric suppression of twist 3 relative to twist 4 we obtain S"^^ <^ 
f^"^^ also numerically. (Generally, a positive feature of our approach is that the parametric 
ordering of quantities is borne out by the numerical values, i.e., that parametrically 
suppressed contributions are also small numerically.) The small value for rf^^ is consistent 
with present experimental data on polarized structure functions. 

The numerical magnitude of the higher-twist matrix elements obtained in our ap- 
proach can be understood in simple terms. The constituent quark picture implied by the 
instanton vacuum means that the typical magnitude of the nucleon matrix elements is 
M'^/Mff times a number of order unity for parametrically large matrix elements (twist 
4), or a number <C 1 for parametrically suppressed matrix elements (twist 3). 

We have found that the fiavor-nonsinglet part of power corrections to polarized struc- 
ture functions, /^g and d}^sy leading in the l/A^'c-expansion. In the unpolarized case 
the situation is opposite, and the flavor-singlet part is leading. This ordering agrees well 
with the results of QCD sum rule calculations of |2^, |3^ 

To summarize, the instanton vacuum implies a definite hierarchy of higher-twist ma- 
trix elements. With increasing accuracy of the measurements of polarized and unpolarized 
structure functions one should be able to test the predictions of this picture more accu- 
rately. 

The methods developed in this paper can be applied to calculate a wide range of other 
matrix elements relevant to DIS. For example, one can compute the twist-3 contribution 
to higher moments of g2, eqs.(|2.8|, p^) , and reconstruct the entire structure function. 
Furthermore, the formalism developed in section ^]2| can be applied to the calculation 
of nucleon matrix elements of 4-fermionic higher-twist operators which appear in 1/ Q^- 
and 1 / Q^-power corrections to unpolarized DIS |^ . 
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(a) 




(b) 



Figure 1: (a) The effective operator, eq.( |3.25|) , in the case of one quark ffavor, Nf = 1. 



(b) Diagrams contributing to the matrix element of the effective operator in quark states. 
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(i) 



(ii) 



( Hi) 



Figure 2: Feynman diagrams contributing to the nucleon matrix element of the effec- 
tive operator, (i) "Valence quark" contribution, eq. (|4.23| ) (ii) "Sea quark" contribution, 
eq. ([4.24|) . (Hi) Interaction contribution, eq.( 4.27|) . 
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